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Abstract 



We analyze the behavior of the ensemble of surface boundaries of 
the critical clusters at T = T c in the 3d Ising model. We find that 
N g (A), the number of surfaces of given genus g and fixed area A, 
behaves as A~ x ^ e~^ A . We show that [i is a constant independent 
of g and x(g) is approximately a linear function of g. The sum of 
Ng(A) over genus scales as a power of A. We also observe that the 
volume of the clusters is proportional to its surface area. We argue 
that this behavior is typical of a branching instability for the surfaces, 
similar to the ones found for non-critical string theories with c > 1. 
We discuss similar results for the ordinary spin clusters of the 3d Ising 
model at the minority percolation point and for 3d bond percolation. 
Finally we check the universality of these critical properties on the 
simple cubic lattice and the body centered cubic lattice. 
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It is a long-standing hope in theoretical physics that it is possible to find 
a formulation of three-dimensional phase transitions dual to the usual order 
parameter field theory. Such a description recasts the dynamics in terms of 
fluctuating surfaces. In the 3d Ising model such a reformulation is indeed 
possible on the lattice [1], and the relevant problem is the existence of a 
continuum limit for the surface theory. Despite much theoretical effort, real 
progress in this direction has been slow and difficult. In order to provide 
some experimental data to perhaps point in the right direction, we embarked 
on a numerical study of surfaces in the 3d Ising model. 

The main obstacle to a continuum surface description is the well-known 
fingering instability of surfaces in embedding space of dimension d > 2 [2]. For 
d < 2 the theory of random surfaces is solved [3], providing a striking confir- 
mation of the scaling predictions of Liouville theory [4]. Another motivation 
for our study was to see if the d > 2 instability plagues Ising surfaces, and to 
look for scaling properties in the distribution of surfaces of fixed genus. We 
will present data that shows nontrivial topology-dependent scaling behavior 
for self-avoiding surfaces in 3d (for previous work, see [5]). 

Contrary to the widely held belief that a phase transition is characterized 
by nested clusters of ordered domains of all possible sizes, it has been shown 
[6] that the distribution of domain boundaries in the 3d Ising model does 
not scale at criticality. Rather, as the critical point is approached from low 
temperatures, islands of flipped spins (which we shall call minority clusters) 
merge into a large percolating cluster at a temperature T p well below the 
critical temperature T c . The dynamics is then dominated not by the entropy 
in the cluster distribution but by the entropy of configurations of the perco- 
lating cluster. Moreover, no local order parameter of the Ising model reflects 
this percolation transition. 

However there exists a cluster representation of the Ising model due to 
Fortuin and Kasteleyn [7] which captures the critical properties of the model, 
in particular the divergence of the correlation length is related to the perco- 
lation of the FK clusters. These clusters are formed by adjoining neighboring 
spins with a probability 1 — e _2/3 if they are equal. The FK representation led 
Swendsen and Wang to propose a Monte Carlo that partially defeats critical 
slowing down [8]. In our simulations we have used the SW algorithm. We 
studied the self-avoiding surfaces bounding both minority spin clusters and 
FK clusters. To explore further the possible realizations of self-avoidance, we 
also simulated the surfaces defined by pure bond percolation [9] . It should be 
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emphasized that the 'bosonic' surfaces defined by Ising domain walls are not 
the same as the 'fermionic' surfaces that arise in the surface reformulation of 
the lattice Ising model; nevertheless we expect that they should capture the 
characteristic features of any critical surface theory of the 3d Ising model. 

We ran a medium sized simulation, using roughly four months of time 
on RISC workstations. We have analyzed Ising configurations on a 64 3 body 
centered cubic (BCC) lattice at a temperature of f3 — .0857 using .3 x 10 6 
iterations. We also collected data on simple cubic lattices of size 32 3 and 
64 3 at = 0.221651 (performing about 6 x 10 6 and .25 x 10 6 iterations 
respectively). Data were also taken at the minority percolation point for the 
Ising model (where we also studied many different lattice sizes, going from 
32 3 to 100 3 ) and for pure bond percolation on the simple cubic lattice. All of 
our error analysis has been done by using jack-knife and binning techniques. 
For more details see ref. [10]. 

On the BCC lattice, we coupled with equal strength both the 6 nearest 
and 8 next-nearest Ising spins so that only three plaquettes of the dual lattice 
meet along a dual link. Since surfaces built this way on the BCC lattice are 
naturally self-avoiding, computing the genus of the dual surface is trivial. 
The number of handles g is obtained through the Euler formula 2 — 2g = 
V—L+P, where V, L, and P are the number of vertices, links, and plaquettes, 
respectively on the dual surface. On the other hand, the genus definition 
on the simple cubic lattice is more problematic, and requires a few choices 
(which we discuss in detail in ref. [10]) to resolve ambiguities where surfaces 
self-touch. 

Overall, for the FK clusters we obtained our best results on the BCC 
lattice. We found a scaling law for the number of clusters of volume V c i, 
N(V c i) ~ V^ T ] the exponent r = 2.22 [11] with a large systematic error 
which could be as high as 3% [10]. We also measured the quantity A d , 
which counts the number of cluster sites on the boundary, and found that 
asymptotically it was proportional to the cluster volume V c i. 

Fig. 1 shows this dependence on the simple cubic lattice. We see that 
for very small volumes, the lattice regularization constrains V c i to equal A c i 
and for intermediate volumes, there is a small deviation from linear scaling 
(as some interior sites begin to appear). The plateau that appears around 
V c i = 3000 indicates the onset of scaling regime where A c i oc V c \ . The growth 
just at the end of the plot is due to the largest cluster, which wraps around the 
lattice and merges with itself to form extra interior points. This plateau is the 
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first indication that the surfaces are not smooth and are unstable towards the 
formation of quasi-one-dimensional objects. The observed proportionality 
of V c i and A d is well-known in the context of pure percolation in 2 and 3 
dimensions [9]. We note that in the well understood case of the 2d Ising model 
(which has a non-pathological continuum behavior) the cluster perimeter 
is not proportional to the area spanned by the cluster. The area-volume 
proportionality, along with the relation [11] to the magnetization exponent 
5 = 1/(t — 2), implies that the scaling behavior of these 'polymers' are related 
to critical properties of the Ising model. 

Turning to the analysis of the topology of the dual surfaces bounding the 
clusters, we consider the distribution N g (A), where g is the genus and A is 
the dual surface area. In Fig. 2, we present our data for genus 5 along with 
a best fit to the functional form 

N g (A) = C g A xi9) e-^ 9)A . (1) 

The fit is superb. One of our main results is that this functional form fits 
our data very well for g > 2 up to about g — 20 where our statistics become 
poor. The 'cosmo logical constant' \x is found to be independent of the genus 
for g > 2 (Fig. 3). The value — 114 ± 3 is proportional to the average 
surface area (in lattice units) per handle. 

In Fig. 4 we plot the exponent x(g) as a function of the genus. Once 
again after a transient region for small genus (g = — 4) we find an almost 
linear behavior in the region g = 5 — 15 with a slope of 1.25 ± 0.1. The 
deviations from linearity are small, and in our observation window they can 
be fitted with an effective exponent correction of order 0.1, or with logarith- 
mic corrections. These make an estimate of the large genus behavior of x(g) 
rather difficult. The region we are observing is still transitory, and we cannot 
exclude that the asymptotic slope of x(g) could equal one at large g. This 
value would be expected if the handles were completely uncorrected as for 
a Poisson distribution P g (A) = ^ 1 (fiA) 9 e~ fJ,A . 

The results presented above nave also been analyzed for the simple cubic 
lattice. The data agree with those of the BCC and point to a good universal 
behavior. The main discrepancies can be traced to the short distance ambi- 
guities which plague the definition of the dual surfaces on the simple cubic 
lattice. The scaling exponents are close for the two lattices. The difference 
among non-universal quantities, e.g. ~ 60 on the simple cubic lattice, can 
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be understood from the ratio of the number of plaquettes of the respective 
Wigner-Seitz cells. 

We will discuss in detail in [10] the analogous features for both the per- 
colation of minority spins of the Ising model and for non-interacting bond 
percolation. The general picture is interestingly the same; indeed, the slope 
dx/dg for bond percolation is compatible with the value for FK clusters. 
However, the slope dx/dg ~ 0.7 ±0.1 in the range g = 3 — 40 for Ising 
minority spin percolation. The same caution as before apply to these slope 
values. Another interesting result is the similarity of the measured area scal- 
ing exponents (r = 2.18±0.05 upon extrapolation to large lattices) as well as 
the linear relation between cluster volume and area at the respective critical 
points (FK clusters at T c and minority spin clusters at T p ). The difference 
of the x(g) estimated for the two interacting theories we are studying is 
interesting. It is possible that in the asymptotic region one will get 1 for 
the slope in both cases, but the large difference of the measured exponents 
(which are, on the contrary, very similar when looking at FK clusters and at 
bond percolation) says at least something about finite size corrections, i.e. 
about the nature of the interaction. Apart from this effect these results indi- 
cate that the cluster distribution scaling is rather insensitive to the average 
cluster density, which is about ten times less for minority spin clusters. At 
least partially we may be observing some universality of different definitions 
of self-avoiding random surfaces. 

We regard the outcome of our topological studies and the behavior V c i ~ 
A d as a strong indication that the cluster boundaries are in a 'branched' 
phase. The topological evidence suggests that the surfaces grow fingers which 
reconnect with a fixed probability per unit area. The cosmological constant 
of the surfaces of fixed genus is nonzero at the critical point. Adjusting the 
temperature away from criticality will only increase /i, as the large surfaces 
are exponentially suppressed (apart from a few surfaces of the size of the 
lattice above the percolation threshold). Therefore there is no relevant pa- 
rameter in the theory that could be tuned to allow large surfaces of low genus. 
One can imagine that some additional parameter (e.g. one that couples to 
the Euler density of the lattice surfaces, which depends on all the spins in 
a fundamental cell) could be fine-tuned to multicriticality (for another ap- 
proach see [12]). Then there would be the possibility to have a scaling theory 
at fixed genus. However, such a tuning will probably not remove the finger- 
ing instability of the surfaces, so we still are faced with the problem that the 
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continuum theory is not a theory of surfaces but of quasi-one-dimensional 
objects. 

The approximate linearity of x with g is strongly reminiscent of the scaling 
behavior of d < 2 random surfaces 1 . But note however that similar scaling 
behavior can be shown in the 'double-scaling' limit of O(N) vector field 
theories [13], so a similar exponent exists for the scaling of random graphical 
networks of fixed topology. 

Finally, in order to get a more accurate description of the geometry of 
the minority clusters, we analyzed the distribution of cluster cross sections 
as a function of their perimeters. Below the critical temperature T c the dis- 
tribution drops off at a scale of the order of a few lattice spacings, providing 
further indication that the surfaces are composed of small, highly intercon- 
nected tubes. Remarkably, around T c (recall that the surfaces have long 
since percolated at T p < T c ) we find very good scaling behavior which must 
be entirely dominated by the cross sections of the large percolated cluster. 
If one were bold enough to view the percolated cluster as describing the 
time evolution of interacting strings, one could speculate that this behavior 
is reminiscent of a Hagedorn-type transition. 

We would like to thank Stephen Shenker for essential discussions which 
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support. After this work was completed, we received an interesting paper 
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ideas related to those presented here. This work was supported in part by 
the Dept. of Energy grants DEFG02-90ER-40560, DEFG02-85ER-40231, 
the Mathematical Disciplines Institute of the Univ. of Chicago, funds from 
Syracuse Univ., by the Centre National de la Recherche Scientifique, by INFN 
and the EC Science grant SC1*0394. 



1 The typical scaling observed for these d < 2 random surfaces is in the free energy, 
describing a single fluctuating surface. Our Ising configurations, however, consist of a 
collection of interacting surfaces. The scalings we observe are for an ensemble of surfaces 
and not the free energy. 
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Figure 1: \n{V c i/A c i) vs. ln(V^) for FK clusters on the L = 64 SC lattice. 



Figure 2: The number of genus 5 surfaces as a function of dual surface area 
A for FK clusters on the L = 64 BCC lattice, with a best fit to the functional 
form given in equation 1. 



Figure 3: The dependence of the cosmological constant fj,(g) on genus for FK 
clusters on the L = 64 BCC lattice. 



Figure 4: The dependence of the exponent x(g) (extracted from the fits to 
equation 1) on genus for FK clusters on the L = 64 BCC lattice. 



10 



0.015 



"i i i i i i i i i i i i i i i i i i i i i i r 



0.01 



0.005 



j j i 







J I I I I I I I I I I I I I I I I I I I I I L 







5 10 15 

genus (L=64 BCC, 3d Ising FK clusters) 



20 



X 



